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A microscopic coupled-channels model for Compton and pion scattering off the nucleon is intro- 
duced which is applicable at the lowest energies (polarizabilities) as well as at GeV energies. To 
introduce the model first the conventional K-matrix approach is discussed to extend this in a follow- 
ing chapter to the "Dressed K-Matrix" model. The latter approach restores causality, or analyticity, 
of the amplitude to a large extent. In particular, crossing symmetry, gauge invariance and unitarity 
are satisfied. The extent of violation of analyticity (causality) is used as an expansion parameter. 



I. INTRODUCTION 



In a K-matrix approach an infinite series of rescattering loops is taken into account with the approximation of 
incorporating only the pole contributions of the loop diagrams. In this approximation lies both the strong and the 
weak sides of this approach. By including only the pole contributions, which correspond to rescattering via physical 
' states, unitarity in obeyed and the infinite series can be expressed as a geometric series and summed as given in 
O . Eq. (§. Such an approach can be formulated in a co-variant approach, where electromagnetic-gauge invariance is 
obeyed through the addition of appropriate contact terms. It can also be shown that crossing symmetry is obeyed. 
. In the structure of the kernel a large number of resonances can be accounted for. The weak point of including only 
' the pole contributions of the loop corrections (i.e. only the imaginary part of the loop integrals) is that the resulting 
amplitude violates analyticity constraints and thus causality. 



<^ ' In Section |^ an extension of the K-matrix model is discussed where, without violating the other symmetries, an 



additional constraint, that of analyticity (or causality), is incorporated approximately. In this approach, called the 
Q . "Dressed K-Matrix Model" , dressed self-energies and form factors are included in the K-matrix. These functions are 
^ ' calculated self-consistently in an iteration procedure where dispersion relations are used at each recursion step to 
. relate real and imaginary parts. 



II. THE K-MATRIX APPROACH 

In K-matrix models the T-matrix is written in the form, 

T = {I - KiS)-^ K , (1) 

where 5 indicates that the intermediate particles have to be taken on the mass shell and all physics is put in the kernel, 
the K-matrix. This amounts to re-sum an infinite series of pole contributions ofloop corrections. It is straightforward 
to check that 5* = 1 -I- 2iT is unitary provided that the kernel K is Hermitian. Since Eq. (|l|) involves integrals only 
over on-shell intermediate particles, it reduces to a set of algebraic equations when one is working in a partial wave 
basis. When both the tt — N and 'y — N channels are open, the coupled-channel K-matrix becomes a 2 x 2 matrix in 
the channel space, i.e. 



K = 



K K 
K K 



(2) 



It should be noted that due to the coupled channels nature of this approach the widths of resonances are generated 
dynamically. 

In traditional K-matrix models the kernel, the K-matrix, is built from tree-level diagrams |, |. In the present 
investigation the type of diagrams included in K^^ are similar to that of ref. [Q except that the A is treated now as 
a genuine spin-3/2 resonance in order to be compatible with the later treatment of the in- medium A resonance. 
This K-matrix is indicated in Fig. (|l|). Most of the (non-strange) resonances below 1.7 GeV have been included. The 
different coupling constants were fitted to reproduce pion scattering, pion photoproduction and Compton scattering 
on the nucleon. A comparable fit to the data as in ref. Q could be obtained. In Fig. (|^) the results for Compton 
scattering are compared to data. 
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FIG. 1: The sum of diagrams included in the K matrix for Compton scattering. A full spectrum of baryon resonances has been 
included. 
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FIG. 2: The calculated cross section for Compton scattering off the proton as a function of angle at fixed photon energy, and 
as a function of photon energy at fixed angle. Data are taken from ref. 



III. BASIC SYMMETRIES 



A realistic scattering amplitude for a particular process should obey certain symmetry relations, such as Unitarity, 
Covariance, Gauge invariance, Crossing and Causality. In the following each of these symmetries will be shortly 
addressed, stating its physical significance. It is also indicated which of these is obeyed by the K-matrix approach 
discussed in the previous section. The comparative success of the K-matrix formalism can be regarded as due to the 
large number of symmetries which are being obeyed. A violation of anyone of these symmetries will directly imply 
some problems in applications. Improvements are thus important. 



A. Unitarity 



The unitarity condition for the scattering matrix S reads 5*^5 — 1. Usually one works with the T- matrix operator 
which can be defined as S = 1 + 2iT, and the unitarity condition is rewritten as 2i{TT'')fi = Tfi — T*j. If the 
T-matrix is symmetric (which is related to time- reversal symmetry), the last formula becomes Im Tfi = X]n-^/"'^in' 
where the sum runs over physical intermediate states. The latter relation is the generalization of the well-known 
optical theorem for the scattering amplitude. Unitarity can only be obeyed in a coupled channel formulation; the 
imaginary part of the amplitude "knows" about the flux that is lost in other channels. 



3 



In the K-matrix formalism the T-matrix is expressed as T ~ i^ix which imphcs that S = ^ is clearly unitary 
provided the the kernel K is Hermitian. This kernel is a matrix, where the different rows and columns correspond to 
different physical outgoing (incoming) channels. The coupled channels nature is thus inherent in such an approach. 
As explained earlier the kernel is usually written as the sum of all possible tree- level diagrams. In a partial- wave basis 
JsT is a matrix of relatively low dimensionality and the inverse, implied in the calculation of the T matrix, can readily 
be calculated. 



B. Covariance 

The scattering amplitude is said to be covariant if it transforms properly under Lorentz transformations. As a 
consequence the description of the reaction observables is independent of the particular reference frame chosen for 
the calculations. It naturally implies that relativistic kinematics is used. 

Since the appropriate four-vector notation and 7-matrix algebra are used throughout our calculation, the condition 
of Lorentz covariance is fulfilled. 



C. Gauge invar iance 

Gauge invariance means that there is certain freedom in the choice of the electromagnetic field, not affecting the 

observables. Its implication is current conservation, VJ = or in four- vector notation, 9^J^ = 0. Using the well 
known correspondence between momenta and derivatives, current conservation can be re expressed as fc^ J*^ = 0. If 
the electromagnetic current obeys this relation it can easily be shown that observables, such as a photo-production 
cross section, arc independent of the particular gauge used for constructing the photon polarization vectors. 

One of the sources for violation of gauge invariance is the form factors used in the vertices. A form factor implies 
that at a certain (short range) scale a particle appears 'fuzzy'. At distances smaller than this scale deviations from 
a point-like structure are important; however in the formulation the dynamics at this short scale is not sufficiently 
accurate. For one thing, the flow of charge at this scale is not properly accounted for, implying violation of charge 
conservation. To correct for this, so-called contact terms are usually included in the iiT-matrix kernel. In the present 
model these contact terms axe constructed using the minimal substitution rules. The corresponding T-matrix, as well 
as the observables, are independent of the photon gauge. 

D. Crossing Symmetry 

Physical consequences of the crossing symmetry are more difficult to explain. It basically means that in a proper 
field-theoretical framework the scattering amplitudes of processes in the so-called crossed channels can be obtained 

from each other by appropriate replacements of kinematics. This assiimcs that the amplitiide can be analytically 
continued from the physical region of one channel to the physical regions of other channels. An example of the crossed 
channels is 7 A/' — »■ ttA'', ttA'' — > and NN "/it. 

Crossing symmetry puts a direct constraint on the amplitude for the case that direct and crossed channels are 
identical, as for example for the processes ttN ttN or jN — > 7^^. In these reactions crossing symmetry leads to 
important symmetry properties of the amplitudes under interchange of s and u variables. Due to the fact that in the 
iiT-matrix formalism the rescattering diagrams which are taken into account have only on-shell intermediate particles, 
it can be shown that the s-u crossing symmetry is obeyed provided that the kernel itself is crossing symmetric. Since 
the latter is the case, crossing symmetry is obeyed. 

E. Analyticity 

Analyticity of the scattering matrix is not really a symmetry. Rather, it requires that the amplitude be an analytic 
function of the energy variable and in particular that it obeys dispersion relations. The physics of analyticity is closely 
related to causality of the amplitude as is illustrated in the following example. 

Assume that a signal is emitted by an antenna at time t = 0. At all subsequent times the signal is given by a 
function F{t) while causality requires that at earlier times there was no signal, F{t < 0) = 0. This signal can be 
Fourier-transformed, f{ui) = °^ dt e^'^^ F (t) to explicitly show its energy or frequency dependence. Note that the 
integration region from t = —00 to t = gives zero contribution due to the causality requirement. This transformation 
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can also be considered for complex values of uj. Since the integration interval runs only over positive values for t the 
Fourier integral exists and is a well behaving function for all complex values of energy uj for which Imiuj) > i.e. it is 
an analytic function. For such a function contour integrals in the complex lo plane can be performed and the function 
obeys the Cauchy theorem which is in this context is usually formulated as a dispersion relation, 



Re/(c.) = - / 



showing that for an analytic function the real and imaginary parts are closely related. For example, if the imaginary 
part of an analytic function is given by the curve on the left-hand side of Fig. the real part of this function is 
given by the right-hand side. 



FIG. 3: An example of the real and imaginary parts of an analytic function which are related through a dispersion relation. 



In the traditionally used K-matrix approaches the analyticity constraint is badly violated. The origin of this is 
explained in the following. 

In a field-theoretical calculation of a scattering amplitude one includes rescattering contributions of intermediate 
particles which are expressed as loop integrals. In Fig. a typical loop contribution to the self energy is shown. 
Ignoring terms in the numerator which are irrelevant for the analyticity properties, the corresponding integral can be 
expressed as 

J(p2) = / d^k -2 2^ = J(p') + tlmJip') (3) 

J [k -- fi + ie\ [[p — ky — 77?/ + it\ 

where the right hand side in this equation and in Fig. (^) expresses the fact this integral has a real and an imaginary 
part, each of which corresponds to some particular physics. The imaginary part of the integral arises from the 
integration region where the denominators vanish, corresponding to four-momenta k where the intermediate particles 
in the loop are on the mass shell -or equivalently- are physical particles with fc^ = and (p—k)^ = w? . Conventionally 
this is indicated by placing a slash through the loop (see Fig. (^) to indicate that the loop can be cut at this place 
since it corresponds to a physical state. The other parts of the integration region contribute to the real part of the 
integral. In the latter case the particles in the loop are off the mass shell. 



FIG. 4: Loop integral contributing to the self energy. 



It can be shown that the K-matrix formulation for the T-matrix corresponds to including only the imaginary (or 
cut-loop) contributions of a certain class of loop diagrams. This guarantees (as was shown before) that unitarity is 
obeyed. Analyticity of the scattering amplitude is however violated due to ignoring the real contributions of these 
loop integrals. As a consequence causality will be violated! 

To (partially) recover analyticity of the scattering amplitude the so-called "Dressed K-matrix approach" has been 
developed. It is described in the following section. 
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IV. THE DRESSED K-MATRIX MODEL 

As discussed in the previous section, the coupled channels K-matrix approach is quite successful in reproducing 
Compton scattering. However it fails in predicting nucleon polarizabilities. The reason is that, in spite of the many 
symmetry properties that are satisfied, analyticity or causality of the amplitude is badly violated. In the "Dressed 
K-matrix" approach the constraint of analyticity is incorporated in an approximate manner without spoiling the other 
symmetries. In fact analyticity is used as a kind of expansion parameter where presently only the leading contributions 
are included. The ingredients of the Dressed K-Matrix Model were described in Refs. |^, ^ and the main results 
were presented in Rcf. The essence of this approach lies in the use of dressed vertices and propagators in the 
kernel K. 

The objective of dressing the vertices and propagators is solely to improve on the analytic properties of the ampli- 
tude. The imaginary parts of the amplitude are generated through the K-matrix formalism (as imposed by unitarity) 
and correspond to cut loop corrections where the intermediate particles are taken on their mass shell. The real parts 
have to follow from applying dispersion relations to the imaginary parts. We incorporate these real parts as real vertex 
and self-energy functions. Investigating this in detail (for a more extensive discussion we refer to |6|) shows that the 
dressing can be formulated in terms of coupled equations, schematically shown in Fig. (^, which generate multiple 
overlapping loop corrections. The coupled nature of the equations is necessary to obey simultaneously unitarity and 
analyticity at the level of vertices and propagators. 

Off-Shell On-shell 




Free N propagator 



FIG. 5: Graphical representation of the equation for the dressed irreducible ttNN vertex, denoted by an open circle, and the 
dressed nucleon propagator, denoted by a solid line. The dashed lines denote pions, the double lines denote As and the zigzag 
and dotted lines are p and a mesons, respectively. The resonance propagators are dressed. The last term in the second equation 
denotes the counter-term contribution to the nucleon propagator, necessary for the renormalization. 

The equations presented in Fig. (||) are solved by iteration where every iteration step proceeds as follows. The 
imaginary - or pole ~ contributions of the loop integrals for both the propagators and the vertices are obtained by 
applying cutting rules. Since the outgoing nucleon and the pion are on-shell, the only kinematically allowed cuts are 
those shown in Fig. (^). The principal- value part of the vertex (i.e. the real parts of the form factors) and self-energy 
functions are calculated at every iteration step by applying dispersion relations to the imaginary parts just calculated, 
where only the physical one-pion-one-nucleon cut on the real axis in the complex p^-plane is considered. These real 
functions are used to calculate the pole contribution for the next iteration step. This procedure is repeated to obtain 
a converged solution. We consider irreducible vertices, which means that the external propagators are not included 
in the dressing of the vertices. 
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Bare ttNN form factors have been introduced in the dressing procedure to regularize the dispersion integrals. The 
bare form factor reflects physics at energy scales beyond those of the included mesons and which has been left out of 
the dressing procedure. One thus expects a large width for this factor, as is indeed the case. 

The dressed nucleon propagator is renormalized (through a wave function renormalization factor Z and a bare 
mass mo) to have a pole with a unit residue at the physical mass. The nucleon self-energy is expressed in terms of 
self-energy functions A{p'^) and B{p^) as Y,iq{p) = A]y{p'^)j/+ Bn{p^) m. 

The procedure of obtaining the ^NN vertex [Q is in principle the same as for the nNN vertex. Contact ^ttNN and 
"fjNN vertices, necessary for gauge invariance of the model, are constructed by minimal substitution in the dressed 
ttNN vertex and nucleon propagator, as was explained in 0|. 

The present procedure restores analyticity at the level of one-particle reducible diagrams in the T-matrix. In 
general, violation due to two- and more-particle reducible diagrams can be regarded as higher order corrections. An 
important exception to this general rule is formed by, for example, diagrams where both photons couple to the same 
intermediate pion in a loop (so-called "handbag" diagrams). This term is exceptional since at the pion threshold the 

1/2 

S-wave contribution is large, due to the non-zero value of the Eq[^_ multipole in pion-photoproduction, leading to a 
sharp near-threshold energy dependence of the related Compton amplitude [p^ . In the K-matrix formalism, the 
imaginary (pole) contribution of this type of diagrams is taken into account. Not including the real part of such a 
large contribution would entail a significant violation of analyticity. To correct for this, the 'yjNN vertex also contains 
the (purely transverse) "cusp" contact term whose construction is described in Section 4 of Ref. [Q. Since, due to 
chiral symmetry, the S-wave pion scattering amplitude vanishes at threshold, the mechanism that gives rise to the 
important "cusp" term in compton scattering does not contribute to 'K'kNN or ttjNN contact terms. The analogous 
to the "cusp" "fjNN term will thus be negligible and have therefore not been considered. 



A. Results 



Results for pion-nucleon scattering and pion-photoproduction obtained in the dressed K-matrix model and in the 
traditional K-matrix approach are of comparable quality. One should, however, expect the two approaches to have 
significant differences for Compton scattering since for this case constraints imposed by analyticity will be most 
important ||l2|, |l3|. 

The effect of the dressing on the amplitude can be seen in Fig. (^, where also the results of dispersion analyses 
are quoted for comparison. Note that the imaginary parts of from calculations B (Bare, corresponding to the 
usual K-matrix approach) and D (Dressed, the full Dressed K-matrix results) are rather similar in the vicinity of 
threshold. 
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FIG. 6: The partial amplitude of Compton scattering on the proton in units lO^^/m^. Solid line: dressed K-matrix, 
D; dotted line: bare K-matrix, B. Also shown are the results of the dispersion analyses of Ref. [^2| (dash-dots) and Ref. p^ ] 
(dashed). 



The polarizabilities characterize response of the nucleon to an externally applied electromagnetic field y_5| 
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They can be defined as coefficients in a low-energy expansion of the cross section or partial amplitudes of Compton 
scattering. Since gauge invariance, unitarity, crossing and CPT symmetries are fulfilled in both models the Thompson 
limit at vanishing photon energy is reproduced. Our results for the electric, magnetic and spin polarizabilitics of 
the proton are given in Table |, where they are compared with the results given in refs. |l^ and with the values 
extracted from recent experiments. The contribution from the the t-channel Tr'^-exchange diagram has been subtracted. 
The effect of the dressing on the polarizabilitics can be seen by comparing the values given in columns D (dressed) 
and B (bare). In particular, the dressing tends to decrease a while increasing (3. Among the spin polarizabilitics, 
7E1 is affected much more than the other 7's. The effect of the additional "cusp" jjNN contact term^, strongly 
influences the electric polarizabilitics rather than the magnetic ones. This is because the "cusp" contact term affects 
primarily the electric partial amplitude /^"^ (corresponding to the the total angular momentum and parity of the 
intermediate state J"^ = l/2~) rather than the magnetic amplitude {J"^ = 1/2^). 

TABLE I: Polarizabilities of the proton. The units are 10~'*/m'^ for a and f3 and 10~''/m'* for the 7's (the anomalous tt" 
contribution is subtracted). The first two columns contain the polarizabilities obtained from the present calculation; D (full, 
dressed) and B (bare K-matrix) . The two columns named xPT contain the polarizabilities calculated in the chiral perturbation 
theory [ p^ , Results of recent dispersion analyses are given in the last column (Ref. for a and /3 and Ref. for the 
7's). 
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Of special interest is to check whether polarizabilities as extracted form the low energy behavior of the amplitude 
are in agreement with the values as extracted from energy weighted sumrules. The derivation of sumrules is based on 
the fact that the amplitudes obey certain symmetries where analyticity is of particular importance. This comparison 
is still in progess and results will be published in a forthcoming paper pO) . Preliminary results indicate that the 
different sum rules are obeyed, with the exception of the sum rule for the spin polarizability 70. This may be due to 
an incomplete dressing of the A-resonance in the present calculational scheme. 

V. SUMMARY AND CONCLUSIONS 

Results are presented for Compton scattering on a free proton as well as on a nucleus. For processes on the proton 
the usual K-matrix model as well as the recently developed dressed K-matrix model are discussed. In the latter 
approach the real self energies and vertex functions are obtained from the imaginary parts using dispersion relations 
imposing self-consistency conditions. It is indicated that such an approach is essential to understand features seen in 
the data, in particular at energies around and below the pion production threshold. 
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